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xhe  objective  of  this  study  was  to  extend  the  SEASAT  altimeter  determi¬ 
nation  oi  the  M2  ocean  tide  from  a  single  point  near  Cobb  seamount  in  the 
northeast  Pacil  Lc  ocean  to  all  the  world's  oceans  on  a  5°  grid.  After 
extending  the  grid  to  32  locations  covering  the  northeast  Pacific  ocean,  it 
bee...  e  apparent  that  the  resulting  tidal  charts  contained  unrealistic  features. 

n  addition,  large  discrepancies  existed  between  the  altimeter  determined  para¬ 
meters  and  tnose  determined  from  bottom  pressure  gauges.  Further  extension 
ot  t  ne  altimeter  solutions  was  there  lore  stopped  and  the  study  efforts  were 
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redirected  Lo  study  the  reason  tor  the  discrepancies. 

It  was  concluded  that  the  uncertainty  of  the  altimeter  solutions  was  :..uci 
higher  than  originally  anticipated,  due  to  the  near-resonance  of  the  SEASAT 
orbit  witli  the  dominant  tidal  components  and  the  relative  shortness  of  the 
bEnbAT  mission.  The  extension  of  the  altimeter  tide  solution  for  M2  over  the 
world's  oceans  requires  data  from  a  future  altimeter  satellite  mission  spanning 
at  least  one  year. 


-  i'riOSMiX 

BACKGROUND 

in  Che  previous  contract  .\UU014-79-C-u409 ,  SEASaT  altimeter  data  was  ana¬ 
lyzed  lor  tidal  information  at  a  location  in  the  northeast  Pacific  ocean  near 
COBB  seamount  (Brown,  [1983J).  The  amplitude  and  Greenwich  phase  angle  of  the 
m2  component  determined  by  harmonic  analysis  of  the  altimeter  sea  height  re¬ 
siduals  were  in  fair  agreement  with  h2  tide  parameters  determined  using  a  bot¬ 
tom  pressure  gauge  on  tobb  seamount  (Larsen  and  Irish,  1975).  The  altimeter 
M2  parameters  are  43  +  8  cm  amplitude  and  243  +  10  degrees  Greenwich  phase 
angle,  while  the  bottom  pressure  gauge  results  an-  81  cm  amplitude  and  242° 
Greenwich  phase.  in  addition,  the  altimeter  solution  seemed  quite  robust. 

The  altimetric  tide  parameters  remained  quite  consistent  as  different  altimeter 
passes  were  deleted  from  the  solution  (see  Figure  1).  The  generally  reduced 
amplitude  of  the  altimeter  solutions  was  thought  to  be  due  to  the  use  of  non- 
optimum  weighting  of  the  different  tidal  components  in  the  solution. 

Encouraged  by  these  results,  in  the  present  contract  we  undertook  to  ex¬ 
tend  tnese  altimetric  tide  solutions  for  M2  over  all  the  world's  oceans,  wher¬ 
ever  sufficient  SEABaT  data  existed  to  permit  the  near  coincident  intersections 
of  5  or  more  north  bound  passes  with  5  or  more  southbound  passes.  In  general, 

this  would  result  in  a  global  grid  of  altimetric  tide  solutions  at  an  interval 
.  o 

of  a  .  At  these  points  the  altimetric  tide  parameters  would  be  compared  to 
those  interpolated  from  bottom  pressure  gauges  or  other  deep  ocean  tide  models. 

aLM'Lafi  mf  aLIIVIIIEB 

By  June  of  19d2,  tile  altimeter  solutions  had  been  extended  to  eleven  loca¬ 
tions  in  the  northeast  Pacil ic  (Brown,  1982)  and  while  5  solutions  agreed  well 
with  nearby  bottom  pressure  gauge  results,  some  systematic  discrepancies  were 
observed  in  the  other  5  solutions.  By  September  of  1982,  the  northeast  Pacific 
grid  had  been  expanded  to  32  solutions,  anti  the  discrepancies  became  more  fre¬ 
quent  and  less  systematic.  These  results  were  described  in  a  paper  presented 
at  UGEAuB  '82  conference,  September  2u  -  22,  1982  in  Washington,  D.G.  As  it 
was  obviously  premature  to  continue  the  global  extension  of  the  solution  in  the 
iace  ol  these  unsolved  problems,  we  decided  to  analyze  the  existing  solutions 
in  uepth .  This  final  report  describes  the  results  of  tiiat  analysis  as  well  as 
a  description  of  Lhe  32  altimetric  tide  solutions  in  the  northeast  Pacific. 
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DERIVATION  OF  SEa.  SURFACE  liEIChTS 

historically ,  the  primary  problem  for  recovery  oi  ocean  tides  from  satellite 
altimetry  is  the  aliasing  of  orbit  error  with  derived  sea  surface  height  (Brown 
and  Hutchinson,  19dl).  Figure  2  illustrates  this  problem. 


Derived  Sea  Surface  height,  S. 

Figure  2 

The  error  £  in  the  knowledge  of  satellite  distance  Rs  from  the  geocenter, 
couples  directly  into  the  computation  lor  se.j  surt.ice  height  S  as  iollows: 

S  =  KS  —  Re  -  h  +  e  (1) 

where  lie  is  radius  to  the  reference  ellipsoid  and  ii  is  the  measured  altitude. 
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The  radial  (height)  component  of  the  satellite  orbit  error  is  modeled  as 
a  low  order  harmonic  function  of  the  orbit  period.  It  is  crucial  for  tile 
separation  ol  tides  from  these  orbit  perturbations  that  the  data  arc  be  longer 
than  a  lull  spatial  wavelength  ol  the  ocean  tide  pattern  ol  interest.  The 
data  arc  length  cnosen  for  orbit  error  removal  was  bOUO  km,  corresponding  to 
a  Liine  span  of  lb  minutes  at  satellite  speed  along  the  subtrack.  The  radial 
orbit  error  e  over  this  time  span  is  modeled  as 


e  =  a  +  BcosW  t  +  Csin  W  t 

+  boos  2W°t  +  Esin2W°t 
o  o 


where  is  the  orbital  frequency  in  radians/sec  and  t<<2it/WQ.  For  the  tidal 

analysis,  we  replace  the  actual  geoid  height  N  by  a  good  mathematical  model 

N  (t),  such  as  GLll- 1011  (Lurch  et  al.,  1978)  at  each  point  of  the  data  arc,  and 
m 

express  the  difference  between  the  derived  sea  height  and  the  geoid  as 


8(t)  -  N(n(t)  =  e(t)  +  d 


where  d  is  unmodeled  geoid  and  temporal  sea  surface  heights  as  well  as  random 
noise.  An  equation  3  is  generated  for  each  altimeter  observation  and  these  are 

is  *  ? k  is 

solved  in  a  least-squares  procedure  for  the  optimal  parameters  A  ,  8  ,  C  ,  D  , 

is 

and  L  .  The  corrected  sea  surface  height, 


s  (l)  =  s ( l 3  -  e  (t)  (a; 

is  ..  least  squares  best  lit  to  the  (JLM-lUb  geoid  profile  along  the  data  arc, 
with  short  wavelength  residuals  due  to  detailed  geoid  anomalies  and  temporal 
ocean  surface  variations. 

DATA  DlSTKiBDTlON 


The  three-day  repeat  orbit  ot  SLaSmI  results  in  a  locally  high  concentration 
of  subtrack  crossings,  oceuring  on  a  regular  grid  spacing  at  about  8  degrees  of 
longitude . 
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Figure  3  displays  108  SEA SAT  arcs  with  crossovers  concentrated  in  32  loc¬ 
ations.  nltimetric  tide  parameters  have  been  determined  for  all  of  these 
locations.  If  we  label  the  northbound  (southeast-northwest)  passes  by  index 
i,  and  the  southbound  passes  by  index  j,  then  orbit-corrected  sea  heights 
interpolated  to  these  crossovers  may  be  written  as 


* 

S  (t .  .)  =  H(t .  .)  +  T(t . .)  +  e. . 

Ji  Ji  Ji 

where  t,  is  the  time  at  which  the  satellite  on  subtrack  1,  overflies  track 

lr.) 

m;  N(t,  )  is  the  actual  geoid  height  at  the  crossover  point;  T(t,  )  is  the 
lm  lm 

height  of  the  ocean  tide  at  the  time  t.  ;  and  e,  is  the  non-tidal  contri- 

Lm  Jm 

bution  to  sea  surface  height  at  this  crossover  at  time  t^m,  such  as  due  to 
eddies,  inverse  barometer  effects,  and  altimeter  noise. 

While  each  of  these  measurements  has  a  unique  location  and  time,  their 
proximity  in  space  and  in  time  (for  a  given  pass),  limits  the  number  of  essen¬ 
tially  different  tide  heights  observed.  For  recovery  of  tidal  parameters,  we 
will  later  assume  that  all  the  different  crossovers  in  a  given  locale  exper¬ 
ience  the  same  tide,  but  for  reducing  random  altimeter  noise  effects,  let  us 
now  consider  that  this  series  of  measurements  represent  statistically  indepen¬ 
dent  measures  of  the  instantaneous  sea  height.  Note  that  t  is  not  equal 
tj^,  and  these  times  may  differ  by  amounts  ranging  from  several  hours  to 
several  months. 

In  general,  it  is  very  difficult  to  extract  tidal  information  from  these 
derived  sea  heights,  because  the  actual  geoid  height  N(t^m)  is  not  accurately 
known.  However,  we  may  apply  the  constraint 

N(tij )  =  (6) 

and  form  difference  equations  which  are  independent  of  knowledge  of  the  actual 
geoid,  as  follows: 
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AS  .  .  =  S*(t..)  -  b*  ( t  .  .  ) 

1J  AJ  J1 

=  T ( t . . )  -  T ( t . .)  +  e'  .  . 

ij  Ji  ij 

where  e '  .  .  =  e .  .  -  e .  . 

ij  J-J  Ji 


(7) 


HARMONIC  ANALYSIS 


We  model  the  oceun  tide  T(t)  for  all  ij  by  a  single  linear  combination  of 
the  5  largest  tidal  harmonics,  assuming  known  period  but  unknown  amplitude  and 
phase.  Implicitly,  we  also  assume  that  ail  the  crossovers  are  collocate.'  each 
observing  tne  same  tide  in  the  local  area.  Substituting  this  expansion  T 
in  equation  7,  we  obtain  , 


As  .  .  =  >  fa,  (sinw.  t  .  .  -sinw.  t .  .  ) 

ij  L—i  L  k  k  ij  k  3 1 


k=l 


+  b,  (cosw,  t .  . -cosw,  t .  .  +  e  '  .  . 

k  k  ij  k  jiJ  ij 

aown  coefficients  rel 
amplitude  and  Greenwich  piiase  angle  4.^  by 


where  the  a.  ,  b.  are  unknown  coefficients  related  to  the  area  mean  tidal 
k  k 


A, 


"  +  lTT 


♦  k  =  Tan"1  ffk\+  K 
bk  / 


where  R  is  the  phase  of  the  kth  component  of  the  equilibrium  tide  at 
K 

Greenwich  at  tile  epoch  time  OOUOhrs  ChY,  Ai>ril  30,  1973  (see  Schureman,  1941). 

even  though  we  are  only  interested  in  the  M2  parameters,  we  parameter ixe 
the  next  largest  components  (kl,  S2,  01,  and  N2)  in  order  to  effect  some  sep¬ 
aration  between  components  which  may  alias  with  M2  in  such  a  short  time  span. 

Thus  we  express  the  difference  in  derived  sea  heights  at  the  subtrack 
crossovers  in  terms  of  differences  in  tidai  phase  angles  of  these  5  components. 
11  the  time  uiflerences  t.  .  -t . .  are  such  that  the  essential  tidal  phase  angle 

YJ  J  L 

differences  are  well-distributed,  then  a  least-squares  solution  of  the  AS,^ 
should  permit  extraction  of  the  tidal  information  from  these  measured  sea 


f <H  OSi'IJX  COil?  Oil  A/  J  Oil 


height  dit f erenees .  For  tidal  components  S2  and  Kl,  the  3  day  repeat  orbit 
does  not  yeild  a  good  distribution  of  tide  phases.  This  is  mitigated  by  the 
use  ot  some  colinear  passes  from  early  in  the  mission  (orbit  number  less  than 
about  1040)  when  the  orbit  track  repeat  period  was  not  so  commensurate  with 
that  of  ocean  tide  components. 

We  seek  a  least  squares  estimation  of  the  ten  tidal  parameters,  a^  and  b^, 
to  minimize  the  sum  of  squares  of  the  quantities 


ASij  “  [T(tij*  ~  T(tji)] 

The  generalized  normal  equation  to  be  solved  is  of  the  form 


(10) 


where  L  is  the  matrix  of  coefficients  such  as 


sinw.  t . ,  -  sinw  t . . 
k  ij  k  ]i 


or 


cosw, t . .  -  cosw. t . . 
k  ij  k  ji 


x  is  the  vector  of  tidal  parameters  a,  ,  b  ,  k  =  1  to  5;  is  the  vector  of 

tC  K 

sea  height  differences  and  W  is  an  arbitrary  weight  matrix  (2K  by  2k). 


We  find  that  a  diagonal  weight  matrix  [W]  is  helpful  in  preventing  ill- 
conditioning,  a  condition  which  causes  the  root-sum-sqt.ured  of  the  solution 
parameters  to  exceed  the  observed  sea  height  differences.  For  a  completely 
unconstrained  solution,  we  set  [W]=0.  If,  on  the  other  hand,  we  desire 
negligible  adjustment  for  certain  parameters,  then  the  associated  weight  matrix 
element  for  the  parameters  to  be  frozen  is  set  to  a  very  large  value.  Since 
we  know  that  the  M2  amplitude  is  generally  larger  than  that  of  the  other 
components,  we  chose  a  constraint  that  influences  the  assignment  of  most  of 
the  tidal  amplitude  to  the  M2  component,  with  lesser  amounts  to  the  other 
components  in  proportion  to  their  relative  amplitudes  in  the  northeast  Pacific 
(Brown,  1983). 
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iliis  type  of  weighting  does  not  force  the  adoption  of  a  particular  value 
for  the  amplitude  of  any  given  tidal  component,  but  merely  guarantees  a  well- 
conditioned  solution.  The  exact  values  of  these  weights  are  not  critical,  and 
have  no  direct  effect  on  the  M2  parameters.  However,  the  phase  angles  of  the 
lesser  components  are  distorted  by  this  form  of  weight  matrix. 


RESULTS 

Altimetric  parameters  for  the  M2  tide  at  all  32  locations  in  the  north¬ 
east  Pacific  are  displayed  in  Figure  4  together  with  those  of  nearby  bottom 
pressure  gauges. 

The  northern-inost  pressure  gauges  are  emplaced  on  seamounts.  From  north 
to  south,  the  Surveyor,  Bowie,  and  Union  seamount  gauges  were  reported  by 
Rapatz  and  iluggett,  [1977],  while  the  Cobb  seamount  gauge  results  are  those  of 
Larsen  and  Irish,  [1975].  The  gauge  west  of  San  Francisco  is  the  OBS-3  sta¬ 
tion  of  Lainont  (Nowroozie,  1972)  as  reported  by  Schwiderski,  [1979].  The 
gauge  furthest  from  the  mainland  is  Josie-II  (Irish,  et  al. ,  1971),  and  the 
souLhurn-most  four  stations  are  Flicki,  Josie,  Kathy,  and  Filloux,  us  reported 
by  Munk,  et  al. ,  [1970].  All  pressure  gauges  are  indicated  by  name,  with  the 
exception  of  station  H,  which  is  a  surface  tide  gauge  reported  by  Schwiderski 
[1979] .  These  bottom  and  surface  tide  gauge  results  are  generally  consistent 
with  the  hypothesis  of  a  counter-clockwise  rotating  amphidromic  system  about 
mid-wav  between  Hawii  and  the  mainland. 

Numerical  values  for  the  altimetric  and  pressure  gauge  parameters  are 
listed  in  Table  I.  The  altimetric  tide  parameters  are  not  strictly  comparable 
with  the  pressure  gauge  results  because  the  altimetric  tide  includes  the  de¬ 
formation  of  the  solid  earth  as  well  as  that  of  the  ocean.  However,  calcu¬ 
lations  by  Parke  and  Hendershott  [1980],  show  that  the  contribution  of  the 
vertical  solid  earth  tide  to  M2  in  this  area  is  less  than  about  5  cm  in 
amplitude  and  essentially  no  change  in  phase.  Thus,  with  this  level  of  pre¬ 
cision,  we  may  directly  compare  bottom  pressure  gauge  and  altimetric  results. 
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Table  1.  Altiraetric  Tide  Parameters 
for  the  Northeast  Pacific 


No.  of 

M2 

Residual  RMS  (cm) 

Locat ion 

Passes 

cm(duK . ) 

Before/After 

H 

_ 

20  (33) 

6 

17 

22(265) 

79/71 

7 

16 

3  (8) 

37/  8 

8 

12 

18(175) 

65/28 

9 

18 

18  (70) 

63/61 

Filloux 

— 

19(107) 

10 

15 

15  (19) 

60/57 

11 

16 

3(135) 

85/33 

12 

13 

6(354) 

25/13 

13 

16 

10(216) 

71/70 

Kathy 

— 

29(128) 

14 

16 

10(142) 

84/83 

Flicki 

— 

43(150) 

15 

15 

9(358) 

43/16 

16 

17 

20(120) 

50/47 

17 

13 

10  (63) 

23/19 

Josie-11 

— 

27(267) 

18 

17 

3(219) 

15/12 

19 

14 

16  (56) 

66/24 

OBS-3 

— 

54(206) 

20 

15 

2(271) 

37/19 

21 

16 

3(339) 

33/32 

22 

13 

3(268) 
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Altimeter  solutions  9,  13,  14,  18,  23,  25,  30,  36,  37,  and  39  yield  fair 
to  good  agreement  with  the  phase  angles  of  nearby  bottom  gauges,  but  the  ampli¬ 
tude  is  generally  underestimated.  Altimeter  solutions  6,  11,  12,  17,  19,  24, 

28,  29,  32,  33,  and  37  seem  to  exhibit  a  phase  angle  lead  of  90  to  180 
degrees  compared  to  nearby  bottom  gauges,  and  again  the  amplitudes  are  under¬ 
estimated.  Altimeter  solutions  7,  10,  15,  20,  21,  22,  26,  27,  and  31  show 
fair  consistency  in  phase  with  the  phase  angles  predicted  by  Schwiderski 
[1979],  shown  as  gray  arrows  in  Figure  4,  while  solutions  8,  and  16  disagree 
by  almost  180  degrees.  Comparing  these  altimeter  solutions  with  M2  parameters 
predicted  by  Schwiderski  [1979]  for  the  same  locations,  the  average  under¬ 
estimate  of  amplitude  is  54  percent  with  an  RMS  phase  difference  of  79  degrees. 
The  altimetric  amplitudes,  contoured  in  Figure  5,  agree  in  general  with  the 
expected  increase  toward  shore  and  away  from  the  amphidrome  center,  as  shown 
in  the  Schwiderski  predicted  amplitudes,  Figure  6. 

Despite  the  general  agreement  in  the  patterns  of  amplitude  and  phase, 
there  remain  many  discrepancies  in  phase,  and  the  region  of  low  amplitudes 
seems  much  too  large,  compared  to  bottom  pressure  gauge  results.  In  particular, 
less  than  10  cm  amplitude  for  M2  for  the  region  of  solutions  28,  29,  32,  33,  and 
36  is  clearly  indicative  of  a  systematic  error  in  the  altimeter  solutions.  The 
systematic  underestimate  of  amplitude  does  not  appear  to  be  due  entirely  to 
excessive  weights  in  the  constraint  matrix,  since  amplitudes  change  only  slight¬ 
ly  when  the  constraining  weights  are  removed. 

DISCUSSION 

The  low  amplitude  problem  appears  to  be  partly  due  to  the  crossover  data 
itself.  A  review  of  the  magnitudes  of  the  crossover  differences  before  recovery 
of  tide  parameters  (see  Figure  7)  shows  a  pattern  very  similar  to  that  of  the 
recovered  M2  amplitudes  (in  Figure  5).  It  is  impossible  to  recover  an  M2  ampli¬ 
tude  of,  say,  80  or  90  cm  as  at  location  32,  from  crossover  differences  whose 
RMS  is  9  cm.  The  RMS  of  crossover  differences,  assuming  these  differences  to  be 
due  entirely  to  a  single  sinusoidal  tide  component,  should  be  indicative  of  the 
peak  amplitude  of  the  tide  component.  Since  the  RMS  of  a  sinuisoid  is  just 
a/VT  times  its  peak  amplitude,  and  the  RMS  of  the  crossover  differences  (the 
square  root  of  the  mean  of  the  squares  of  the  difference  between  two  sinusoidal 
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signals  of  the  same  amplitude  but  different  phases)  is  yr  times  the  RMS  of  a 
single  sinusoidal  process,  the  crossover  RMS  should  ideally  be  equivalent  to 
the  peak  amplitude.  Of  course  this  assumes  no  noise  and  sufficient  different 
samples  of  phase  for  accurate  representation  of  the  sinusoid.  Using  the  rule- 
of-thumb  that  the  peak  value  of  a  noisy  set  of  data  is  generally  4  times  lar¬ 
ger  in  magnitude  than  the  RMS  value  of  the  data  set,  we  assume  that  our  sparse 
sampling  of  phase  differences  permits  a  maximum  possible  M2  amplitude  of  about 
4  times  the  RMS  of  the  crossover  differences.  This  further  assumes  that  only 
M2  contributes  to  the  crossover  differences.  By  this  rule,  the  observations 
for  location  32  seem  to  preclude  recovery  of  the  M2  amplitude  of  77  cm  pre¬ 
dicted  for  this  location  (Schwiderski,  [1979]).  However,  this  may  be  due 
to  an  unfortunate  distribution  of  observations. 

What  then  is  the  reason  for  such  small  recovered  amplitudes  at  other  loc¬ 
ations?  A  detailed  analysis  of  the  crossover  differences  and  their  average 
M2  phase  differences  reveals  a  randomness  in  phase.  For  example.  Figures  8  and 
9  show  the  residuals  for  locations  18  and  39,  together  with  the  sinusoidal 
signals  which  result  from  the  altimetric  parameters  and  from  the  M2  parameters 
predicted  by  Schwiderski  [1979],  For  location  18,  Schwiderski  predicts  an 
amplitude  of  39  cm  and  a  Greenwich  phase  of  249  degrees.  While  the  crossover 
differences  at  location  18  are  of  sufficient  magnitude  to  support  an  M2  ampli¬ 
tude  of  this  size,  their  phases  show  a  randomness  that  the  harmonic  analysis 
interprets  as  near-zero  amplitude  at  the  M2  period.  Even  when  the  agreement 
between  the  altimeter  parameters  and  Schwiderski ' s  predicted  parameters  is 
good,  a3  at  location  39  shown  in  Figure  9,  this  result  is  largely  due  to  good 
fortune.  The  residuals  show  a  great  deal  of  randomness  at  the  M2  period,  and 
the  recovered  M2  parameters  are  supported  by  only  a  few  of  the  crossover 
differences.  The  same  situation  holds  for  all  32  altimetric  tide  solutions. 

Since  the  only  error  source  that  can  be  characterized  as  random  in  this 
process  is  that  of  the  altimeter  itself,  and  since  the  magnitude  of  this  error 
is  only  fr-15  cm  RMS  (Townsend,  [1980]),  the  observed  randomness  at  the  M2 
period  must  be  due  to  aliasing  with  the  other  harmonic  processes  affecting  the 
derived  sea  height.  These  are  the  orbital  error,  with  a  period  of  3/43  days, 
and  the  other  tidal  components,  with  periods  very  close  to  1  day  and  1/2  days. 


Location  18 


Sea  Height  Residual  (cm) 
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If  the  contribution  of  these  sources  to  the  crossover  difference  could  be 
removed  exactly,  then  the  residuals  shown  in  Figures  8  and  9  would  presumably 
appear  much  less  random.  Unfortunately,  we  don't  have  enough  different  phase 
samples  to  allow  reliable  separation  of  these  competing  harmonic  processes. 

For  a  location  with  15  passes,  there  are  essentially  only  14  different  phase 
differences  at  any  given  frequency.  To  separate  the  5  competing  tidal  com¬ 
ponent  periods,  we  should  have  50  to  100  well-distributed  phase  differences  for 
each  period.  This  may  well  require  several  hundred  passes  of  data.  Since  the 
SEASAT  mission  lasted  only  100  days,  we  will  need  at  least  a  year  or  two  years 
worth  of  altimeter  data  from  a  future  altimeter  satellite  mission  to  complete 
the  global  altimetric  tide  solution. 

It  is  possible  that  the  randomness  is  partly  due  to  incomplete  removal  of 
orbit  error.  In  particular,  there  is  a  possibility  of  systematic  bias  differ¬ 
ences  between  North-bound  and  South-bound  passes  which  could  alias  with  t lie 
tidal  signals.  This  possibility  was  examined  and,  while  North-South  bias  dif¬ 
ferences  were  found  ranging  from  0  to  81  cm  with  an  average  of  39  cm,  removal 
of  these  biases  change  of  the  M2  parameters  (less  than  10  cm  in  amplitude  and 
10°  in  phase.  Apparently  the  N-S  bias  does  not  alias  strongly  with  the  M2  period. 

CONCLUSIONS  ANO  RECOMMENDATIONS 

The  altimetric  tide  results  for  M2  exhibit  fair  agreement  with  the  bottom 
pressure  gauge  results  at  Cobb  seamount,  and  at  10  of  21  locations  near  bottom 
pressure  gauges.  However,  phase  recover  at  13  of  the  32  locations  is  in  error 
by  greater  than  90  degrees,  and  the  amplitudes  of  the  M2  tide  are  grossly  under¬ 
estimated  for  all  but  2  or  3  locations.  These  poor  results  are  thought  to  be 
due  to  aliasing  between  the  M2  tide  and  other  tide  components.  This  aliasing 
might  be  corrected  by  additional  observations  at  different  phase  angles,  but 
there  are  no  more  SEASAT  data. 

The  consistency  of  the  altimeter  solutions  shows  promise  for  future  im¬ 
proved  remote  measurement  of  ocean  tides.  There  is  still  a  great  deal  of 
room  for  improvement  in  the  quantity  and  quality  of  the  altimeter  data,  the 
choice  of  orbit  repeat  periods,  and  the  techniques  for  orbit  correction  and 
harmonic  analysis.  It  is  heartening  that  such  good  results  have  been  obtained 
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with  the  rudimentary  harmonic  analysis  used  hero.  Perhaps  altimeter  data 
from  a  future  satellite  mission  such  as  CEOSAT  will  he  sufficiently  well  dis¬ 
tributed  in  time  and  of  sufficient  quantity  and  duration  to  prove  the  inde¬ 
pendent  recovery  of  ocean  tides  without  exact  knowledge  of  the  satellite  orhit 
or  local  ocean  tide  parameters. 
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